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ABSTRACT 

We present an analysis of broad emission lines observed in moderate-luminosity active galactic nuclei 
(AGNs), typical of those found in X-ray surveys of deep fields, with the aim to test the validity of single- 
epoch virial black hole mass estimates. We have acquired near-infrared (NIR; 1.0-1.8 /im) spectra of AGNs 
up to z ~ 1.8 in the COSMOS and Extended Chandra Deep Field-South Survey, with the Fiber Multi-Object 
Spectrograph (FMOS) mounted on the Subaru Telescope. These low-resolution (R ~ 600) NIR spectra provide 
a significant detection of the broad Ha emission line that has been shown to be a reliable probe of black hole 
mass at low redshift. Our sample has existing optical spectroscopy (through programs such as zCOSMOS) 
which provides a detection of Mg II, a broad emission line typically used for black hole mass estimation at 
Z > 1, We carry out a spectral-line fitting procedure using both Ha and Mg II to determine the virial velocity 
of gas in the broad line region, the monochromatic continuum luminosity at 3000 A, and the total Ha line 
luminosity. With a sample of 43 AGNs spanning a range of two decades in luminosity (i.e., XL\ ~ lO 44 " 46 ergs 
s -1 ), we find a tight correlation between the rest-frame ultraviolet and emission-line luminosity with a distri- 
bution characterized by (log(AL3ooo/£Ha)) = 1-52 and a dispersion a = 0.16. There is also a close one-to-one 
relationship between the FWHM of Ha and of Mg II up to 10000 km s" 1 with a dispersion of 0. 14 in the dis- 
tribution of the logarithm of their ratios. Both of these then lead to there being very good agreement between 
Ha- and Mg Il-based masses over a wide range in black hole mass (i.e., M BH ~ 10 7_9 M Q ). We do find a small 
offset in Mg Il-based masses, relative to those based on Ha, of +0.17 dex and a dispersion a = 0.32. In general, 
these results demonstrate that local scaling relations, using Mg II or Ha, are applicable for AGN at moderate 
luminosities and up to z ~ 2. 

Subject headings: black hole physics — galaxies: active — quasars: emission lines 



1. INTRODUCTION 

While the study of galaxy evolution has made important 
strides in recent years by being able to weigh individual galax- 
ies (i.e., determine a mass), the field of quasar research is 
grappling with the issue of how to measure accurately the 
masses of a supermassive black holes (SMBHs) for the dis- 
tant quasar population. Here the challenge is greater due to 
the fact that the sphere of influence of a SMBH can only be 
resolved for a limited sample of nearby galaxies whereas the 
dynamical mass of a galaxy can easily be measured due to its 
large spatial extent. A significant leap forward in our abil- 
ity to both accurately and efficiently measure the masses of 
SMBHs, Mbh, at all redshifts will likely lead to new insights 
on questions such as how are black holes fueled, what it is the 



Electronic address: kenta@astro.snu.ac.kr 

1 Department of Physics and Astronomy, Seoul National University, 599 
Gwanak-ro, Gwanak-gu, Seoul 151-742, Republic of Korea 

2 Research Center for Space and Cosmic Evolution, Ehime University, 2-5 
Bunkyo-cho, Matsuyama 790-8577, Japan 

3 Kavli Institute for the Physics and Mathematics of the Universe (Kavli 
IPMU), The University of Tokyo, 5-1-5 Kashiwanoha, Kashiwa 277-8583, 
Japan 

4 The Hakubi Center for Advanced Research, Kyoto University, Yoshida- 
Ushinomiya-cho, Sakyo-ku, Kyoto 606-8302, Japan 

5 Department of Astronomy, Kyoto University, Kitashirakawa-Oiwake- 
cho, Sakyo-ku, Kyoto 606-8502, Japan 

6 National Optical Astronomy Observatory, 950 North Cherry Aveue, Tuc- 
son, AZ 85719, USA 

7 Institute for Astronomy, University of Hawaii, 2680 Woodlawn Drive, 
Honolulu, HI 96822, USA 

8 Departamento de Astronorma, Universidad de Conception, Casilla 160- 
C, Conception, Chile 

9 Astronomical Institute, Tohoku University, 6-3 Aramaki, Aoba-ku, 
Sendai 980-8578, Japan 

10 Max-Planck-Institut fur extraterrestrische Physik, Giessenbachstrasse 1, 
85748 Garching bei Munchen, Germany 

11 INAF - Osservatorio Astronomico di Roma, Via di Frascati 33, 00040 
Monte Porzio Catone, Italy 

12 Department of Astronomy and Astrophysics, The Pennsylvania State 
University, University Park, PA 16802, USA 

13 Max-Planck-Institut fiir extraterrestrische Physik, Postfach 1312, 84571 
Garching bei Munchen, Germany 

14 Dipartimento di Fisica e Astronomia, Universita degli Studi di Bologna, 
viale Berti Pichat 6/2, 40127 Bologna, Italy 

15 INAF - Osservatorio Astronomico di Bologna, Via Ranzani 1, 40127 
Bologna, Italy 

16 NASA/JPL Spitzer Science center, California Institute of Technology, 



1200 East California Boulevard, Pasadena, CA 91 125, USA 

17 Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cam- 
bridge, MA 02138, USA 

18 The Observatories of the Carnegie Institution for Science, 813 Santa 
Barbara Street, Pasadena, CA 91101, USA 

19 Institute for Astronomy, ETH Zurich, Wolfgang-Pauli-Strasse 27, 8093 
Zurich, Switzerland 

20 European Southern Observatory, Karl-Schwarzschild-Strasse 2, 85748 
Garching bei Munchen, Germany 

21 University of California Observatories/Lick Observatory and Depart- 
ment of Astronomy and Astrophysics, University of California, Santa Cruz, 
CA 95064, USA 

22 Department of Physics, Durham University, South Road, Durham DH1 
3LE, UK 

23 Institute of Astronomy and Astrophysics, Academia Sinica, Taipei 
10617, Taiwan 



2 



MATSUOKA ET AL. 



connection with its host galaxy, and how do SMBHs evolve 
within a cosmological framework. 

Spectroscopy enables us to probe the kinematics of ionized 
gas within the vicinity of a SMBH in distant active galac- 
tic nuclei (AGNs) and luminous quasars to infer their black 
hole masses. Traditionally, emission lines (e.g., C IV, Mg II, 
H/3, and Ha) detected in the optical and velocity-broadened 
between 2000 - 20000 km s -1 are used to probe the gravi- 
tational potential well of a SMBH. This lower limit on the 
velocity width has been set somewhat arbitrarily since there 
exists a well-known population o f both type 1 AGNs having 
narrower line widths (i.e., NLS 1 : lOsterbrock & Pogge lll987j) 
and those w ith intermediate-mass black holes dGreene & Hoi 
2004 l2007t) . 

There are methods to determine the luminosity-weighted 
radial distance between the broad-line region (BLR) and cen- 
tral source, Rblr, for AG N (z S, 0-4) through reverberation- 
mapping campaig ns (e.g jBahcall et al.ll 1 972h iPetersonll 1 993h 
iBentz et al. 2009) based on Balmer lines. Even with the com- 
plex nature of the BLR, this characteristic radius is tightly 
correlated with its luminosity (Kas pi et al.l l2000, 2005j), thus 
providing a means to infer such a distance to the BLR in 
large quasar samples based solely on luminosity. Then cou- 
pled with velocity information provides a viral mass esti- 
mate based on a single-epoch spectrum. Such techniques 
have been applied to large quasar samples most notably 
the Sl oan Digital Sky Survey (SPSS). A number of studies 
(e.g.. IVestergaard & Osmerl l2009l ISteinhardt & Elvisl 1201 Ol 
IShen et alJ201 ll) clearly demonstrate that such samples effec- 
tively probe SMBHs above M B h > 1O 9 M at z ~ 2 (an epoch 
of maximal black hole activity) due to the wide area coverage 
and shallow depth. It is important to keep in mind that these 
black hole masses are based on calibrations using lower lu- 
minosity AGNs at low redshift; their application to luminous 
quasars at hig h redshift is not w ell solidified with reverbera- 
tion mapping dKaspi et alJl2007h . 

Deep surveys, such as COSMOS, GOODS and AEGIS, 
are effective probes of black hole accretion at lower masses. 
Given that the b lack hole mass fun ction is steeply declining at 
logM BH > 9 dShen & Kellvll2012l) , studies of the global pop- 
ulation with SDSS are suscept ible to large uncerta inties when 
extrapolated to lower masses (iKellv & Shenl2012l) . While no- 
ble attempts have been made to characterize the low-mass 
end QogM RH < 9: iTamur a et al.l l2006t iHopkins et al.l 12007b 
Merloni & Heinz 2008), these studies have been based on lu- 
minosity and do not consider virial velocities. These deep 
survey fields have considerable X-ray coverage that can be 
utilized to select AGNs that mitigate biases incurred by host 
galaxy dilution and obscuration. Such selection then has the 
potential to effectively probe the lower luminosity AGN popu- 
lation that may be powered by lower mass black holes or those 
accreting at sub-Eddington rates. Followup optical spectro- 
scopic observations are enabling single-epoch virial black 
hole mass estimates (down to 10 7 M Q ) based on the proper- 
ties of their broad emission fines and c ontinuum luminosity 
dTrump et al .1120091: iMerloni et al.ll2010l) . 

Mass estimates for these higher redshift AGNs, which actu- 
ally constitute the majority of the population in deep surveys 
such as COSMOS, rely on Mg II or C IV since the H/3 line 
(used to calibrate recipes based on local samples) is no longer 
available in the optical window at z > 0.8. While the assump- 
tion that the Mg II line i s produced from the same physical 
region as H/3 may hold dMcLure & Jarvisl 120021: IShen et all 
2008), there are studies that indicate that the physics of 



the broad-line region is not so simple dWang et al.l 120091: 
Trakh tenbrot & Netzerll2012l). especially for the most lumi- 
nous quasars (pyiarziani et aD, 120091; [Steinhardt & Silvermanl 
1201 UlMarziani et al.ll20L3T that can have significant outflows 
possibly in response to a more intense radiation field. Further- 
more, it is non-trivial to disentangle the broad Mg II line from 
Fe emission that sits at its base, especially for lower mass 
black holes and maybe even those at the high mass end. 

We present first results of a near-infrared spectroscopic 
survey of broad-line AGNs (BLAGNs) primarily in the 
COSMOS field using the Fiber Multi-Object Spectrograph 
(FMOS) mounted on the Subaru Telescope. With FMOS, 
we now have the capability to simultaneously acquire near- 
infrared spectra of ~ 200 targets over a field of view of 0.19 
square degrees. In this study, we report on the comparison 
between the Ha emission line profile, detected in the near- 
infrared, with that of Mg II present in previously available 
optical spectra. We aim to establish how effective recipes 
(established locally) are to measure single-epoch black hole 
masses out to z ~ 1.8 and for BLAGNs of lower luminosity 
(i.e., lower black hole mass), as compared to those found in 
the SDSS. Our sample is supplemented with AGNs from the 
Extended Chandra Deep Field-South (ECDF-S) Survey that 
reach fainter depths, in X-rays, than the C/ianiira/COSMOS 
survey and improve our statistics at lower black hole masses. 
In Section [2 we fully describe the FMOS observations in- 
cluding target selection, data reduction, and success rate with 
respect to detecting the Ha emission line. We describe our 
method for fitting broad emission lines in Section [3] Our re- 
sults are described in Section |4] Throughout this work, we 
assume H = 70 km s" 1 Mpc" 1 , f2 A = 0.7, S1 M = 0.3, and AB 
magnitudes. 

2. NEAR-INFRARED SPECTROSCOPY WITH 
SUBARU/FMOS 

The capability of FMOS dKimura et alj|2010h to simulta- 
neously acquire near-infrared spectra for a large number of 
objects o ver a wide field offers great potential fo r studies of 
galaxies (lYabe et al Lll20H iRoseboom etaT]l20 1 2l) and AGNs 
dNobutaetalJl2012l) at high redshift. Over a circular region of 
30' in diameter, it is possible to place up to 400 fibers, each 
with a 1 ."2 aperture, across the field. To detect emission lines 
in AGNs over a wide range in redshift, we elect to use the low- 
resolution mode that effectively covers two wavelength inter- 
vals of 1 .05 - 1 .34 fim (J band) and 1 .43 - 1 .77 //m (H band) 
simultaneously. The spectral resolution is A/ A A w 600, thus 
a velocity resolution FWHM ~ 500 km s" 1 at A = 1.5 /xm, 
suitable for the study of broad emission lines of AGN. Un- 
fortunately, this mode requires an additional optical element 
(i.e., VPH grating) that reduces the total throughput to ~ 4% 
at 1.3 /xm and impacts the limiting depth reachable in a few 
hours of integration time. 

Accurate removal of the bright sky background when ob- 
serving from the ground with near-infrared detectors remai ns 
the primary ch allenge. With FMOS, an OH-airglow su p- 
pression filter (Maihara et al. 1993; Iwamur o et al.l l200lb is 
built into the system that significantly reduces the intensity 
of strong atmospheric emission lines that usually plague the 
J band and H band. Equally important, there is the capability 
(cross-beam-switching; CBS hereafter) to dither targets be- 
tween fiber pairs effectively measuring the sky background 
spatially close to individual objects and through the same 
fibers as the science targets. In this mode, two sequential ob- 
servations are taken by offsetting the telescope by 60" while 
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keeping the target within one of the two fibers. The trade off 
is that only 200 fi bers are available in this CBS mode. We 
refer the reader to Kimu raet al.l (1201 Ol) for full details of the 
instrument and its performance. 

Below, we briefly describe our observing program using 
FMOS including the selection of type 1 AGNs, data reduc- 
tion, and success rate with respect to the detection of broad 
emission lines. Complete details of our program will be pre- 
sented in Silverman et al. (in prep.) along with the full catalog 
of emission-line properties of broad-line AGNs in the COS- 
MOS and ECDF-S. 

2.1. Target selection 

Our primary selection of AGNs is base d on their hav- 
ing X-ray emissio n detected by Chandra (lElvis et al.l l2009t 
ICivano et al.ll2.012b within the central square degree of COS- 
MOS (hereafter C-COSMOS). The high surface density of 
AGNs (~ 2000 per deg 2 ) at the limiting depths (/0.5-2.0 kev > 
2 x 10" 16 ergs cm -2 s _1 ) of C-COSMOS field ensures that we 
make adequate use of the multiplex capabilities of FMOS. In 
addition, two FMOS pointings were observed further out from 
the center of the COSMOS field thus we relied upon the cata- 
log of opti cal and near-inf rared counterparts to XMM-Newton 
sources dBrusa et al.ll2010t) for AGN selection. 

We further require t hat optical spectroscopy (Trum pet al.l 
20071; iLillv et al.l2009t) is available for each source that yields 
a reliable redshift and detection of at least a single broad 
(FWHM > 2000 km s" 1 ) emission line, namely Mg II in many 
cases. We then specifically targeted those with spectroscopic 
redshifts that allows us to detect either H/3 (1.2 < z < 1 .7 and 
1.9<z<2.6), Ha (0.6 < z < 1.0 and 1.2 < z < 1.7), or Mg II 
(2.8 < z < 3.8 and 4.1 < z < 5.3) in the observed FMOS spec- 
tral windows in low-resolution mode, i.e., 1.05-1.34 /im and 
1 .43 - 1 .77 /jm. Fibers are assigned to BLAGNs (for which 
we can detect emission lines of interest) with a limiting mag- 
nitude of 7ab = 23. Those at Jab < 21 .5 are given higher pri- 
ority to ensure that this sample (of lower density on the sky) 
is well represented in the final catalog; this also effectively 
improves upon our success rate of detection both continuum 
and line emission. Due to the sensitivity of FMOS and the 
low number density of AGNs at z > 3, our sample has very 
few detections of Mg II in the near infrared. 

In addition, we have acquired F MOS observations o f X- 
ray selected AGNs in the ECDF-S dLehmer et al.ll2005l) sur- 
vey with the inclusio n of those that are only detected in the 
deeper 2-4 Ms data dLuo et alJl2010t IXue etal][2011h in the 
central region that covers GOODS. This deeper X-ray field 
offers the potential to extend the dynamic range of our study 
in terms of black hole mass and Eddington ratio. We specif- 
ically select AGNs, as mentioned above, based on their op- 
tical pr operties determined t hrough deep spectrosc opic cam- 
paigns (ISzokoly et al.ll2004t ISilverman etaDl2010l) . As with 
the COSMOS sample, we place higher priority on the brighter 
AGNs (Rab < 22) while targeting the fainter cases with lower 
priority. 

2.2. Observations 

We have acquired near-infrared spectra with Subaru/FMOS 
of broad-line AGNs from the COSMOS and ECDF-S sur- 
veys. The majority of the data was obtained during open use 
time through NAOJ over three nights in December 2010 (ID; 
S10B-108) and two nights in December 201 1 (ID; S11B-098). 
Additional targets were observed during other programs be- 
ing carried out in the COSMOS field through the University 
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FIG. 1 .— Selection and completeness of the COSMOS AGNs with FMOS 
spectra. Each panel shows the distribution of the observed sample (open 
histogram) and those that yield black hole mass estimates based on the Ha 
emission line (filled black histogram) for (a) X-ray 0.5 - 10.0 keV flux and 
(b) J-band magnitude. 

of Hawaii in S10B-S11A. Weather conditions were accept- 
able although clouds, mainly cirrus, reduced our observing 
efficiency. The typical seeing was ~ 1" with considerable 
variation across the nights. 

We elected to use the CBS mode while taking two sequen- 
tial exposures each of 15 minutes for each position (namely 
A and B positions hereafter). These pairs of exposure were 
repeated multiple times in order to reach an effective total in- 
tegration time of 2.0-3.5 hours on-source. Some time is lost 
to refocusing and repositioning fibers at regular intervals dur- 
ing the full observation. In the early data, only one spectro- 
graph (IRS1) was available thus ^100 fibers were available 
for science targets. 

2.3. Data Reduction 

We use the publicly available sof tware FIBRE-pac ( FMOS 
Image-Based REduction package; Iwamuro et al. 2012j). The 
reduction routines are based on IRAF tasks although several 
steps are processed by additional tools written by the FMOS 
instrument team. Since our observation are carried out us- 
ing an ABAB nodding pattern in CBS mode of the telescope, 
an effective sky subtraction (A-B) can be performed using the 
two different sky images: A„- B„_i and A„- B„+i taken before 
and after the n-th exposure. After the initial background sub- 
traction, a cross talk signal is removed by subtracting 0.15% 
for IRS 1 and 1 % for IRS2 from each quadrant. The differ- 
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ence in the bias between the quadrants is corrected to make 
the average over each quadrant equal. We further apply a flat 
field correction using a dome lamp exposure. Bad pixels are 
masked throughout the reduction procedure. Additional steps 
include the distortion correction and the removal of residual 
airglow lines. This procedure is carried out for both positions 
A and B. Individual frames are combined into an averaged 
image and an associated noise image. Finally, the wavelength 
calibration is carried out based on a reference image of a Th- 
Ar emission spectrum. Individual one-dimensional science 
and error spectra are extracted both to be used for the fitting 
of emission line profiles. 

We perform a first attempt at flux calibration by using spec- 
tra of bright stars, usually 1-2 per spectrograph. A single 
stellar spectrum for each spectrograph is chosen to apply a 
correction based on the spectroscopic magnitude and photom- 
etry from 2MASS. An improvement of the absolute flux level 
is required to account for aperture effects. Therefore, we scale 
the flux of each FMOS spectra of our AGNs to match the deep 
infrared photometry using total magnitudes av ailable from the 
UltraVISTA survey (M cCracken et al.l 1201 2|) available over 
the COSMOS field. While scale factors can reach as high 
as ~ A, the median value is 1 .64. 

2.4. Sample characteristics and completeness 

We have observed over 100 type 1 AGNs in the combined 
COSMOS and ECDF-S fields to date. In Figure Q] we show 
the X-ray flux and NIR magnitude distribution of the 108 
AGNs (originally identified as Chandra X-ray sources) in the 
COSMOS field that have 0.6 < z < 1.8, a redshift interval 
where we are capable of detecting Ha in the FMOS spec- 
troscopic window. The distributions are shown for both the 
observed objects and the 56 having a significant detection of 
a broad Ha emission line (at the expected wavelength) that 
subsequently yielded a black hole mass estimate. 

Prior to our observations, we had no empirical assessment 
of the performance of FMOS thus AGN were targeted to faint 
infrared magnitudes, now understood not to be feasible using 
the low-resolution mode for the faintest objects (7ab ~ 23). 
As shown in Figure \V[b), we have a reasonable level of suc- 
cess (71%) with the detection of broad emission lines at 
brighter magnitudes (Jab < 21.5). Unfortunately, the success 
rate is significantly lower at fainter magnitudes thus dropping 
to 50% for the entire sample with J AB < 23. It is worth high- 
lighting that a survey depth of 7ab < 21 .5 is about two magni- 
tudes fainter than curre nt near-infrared sp ectroscopic observa- 
tions of SDSS quasars (Shen & Liu 2012) at similar redshifts. 
In Figure|2] we demonstrate this by plotting the bolometric lu- 
minosity (based on L3000 and a bolometric correction of 5.15; 
IShen etalJfeOllh of our AGN compared to those from SDSS 
s urveys. 

iGreene & Hoi d2005l) describe in detail the benefits of using 
Ha for black hole mass measurements. In particular, the Ha 
emission line is stronger (~ 3 x) thus more easily detected as 
compared to H/3. This is clearly evident from our observa- 
tions. We successfully detect a broad Ha line in ~ 50% of the 
cases (as mentioned above) while the detection of H/3 is very 
low (17%). Even so, we do detect H/3 in a fair numberof cases 
up to z ~ 2.6 that will be presented in the full emission-line 
catalog (Silverman et al. in prep.). Any future FMOS cam- 
paign designed to detect the H/3 emission line (both broad 
and narrow) should increase the exposure time significantly 
and/or use the high-resolution mode that has roughly three 
times higher throughput in the H band as compared to low- 
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FIG. 2. — Bolometric luminosity as a function of redshift for AGNs in 
COSMOS (red filled circles), ECDF-S (red open circles), and SDSS (grey 
and black symbols). 

resolution mode (see Figure 19 of lKimura et aHl20Tol) . 

3. SPECTRAL LINE FITTING 

We measure the properties of broad emission lines present 
in optical and near-infrared spectra to derive single-epoch 
black hole mass estimates. For the emission lines of inter- 
est here (i.e., Mg II and Ha), we specifically measure the 
full width at half maximum (FHWM), total luminosity of the 
emission line in the case of Ha, and continuum luminosity at 
3000 A. Due to the moderate luminosities of the AGN sam- 
ple, there can be a non-negligible host galaxy contribution that 
impacts the estimate of the AGN continuum at redder wave- 
lengths; therefore, we chose to use the Ha line flux rather than 
the continuum luminosity. Fortunately, the multi-wavelength 
photometry of the COSMOS, including the HST imaging, en- 
ables us to determine the level of such contamination that will 
be fully assessed in a future study. Emission lines are fit using 
a procedure as outlined below that enables us to characterize 
the line shape for even those that have a considerable level of 
noise. 

Our fitting procedure of the continuum and line emission 
utilizes MPFITFUN, a Levenberg-Marquardt least squares 
minimization algorithm as available within the IDL environ- 
ment. Even though this routine has well-known computa- 
tional issues, this algorithm is widely used due to its ease of 
use and fast execution time. The routine returns best-fit pa- 
rameters and their errors as well as measure of the goodness 
of the overall fit. We further describe the individual compo- 
nents required to successfully extract a parameterization of 
the broad component used in determining virial masses. It is 
worth recognizing that each line has its own advantages and 
disadvantages that need to be considered carefully especially 
when fitting data of moderate signal-to-noise ratio (S/N). A 
final inspection of each fit by eye is performed to remove ob- 
vious cases where a broad component is not adequately deter- 
mined almost exclusively due to spectra having low S/N. 

3.1. Ha 

We perform a fit to the Ha emission line (if detected within 
the FMOS spectral window) in order to measure line width 
and integrated emission-line luminosity. Based on the spec- 
troscopic redshift as determined from optical spectroscopy, 
we select the spectrum at rest wavelengths centered on the 
emission line and spanning a range that enables an accurate 
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FIG. 3. — Examples of our fitting routine for the Ha line detected in COS- 
MOS AGN. In the upper portion of each panel, the observed spectrum is 
shown in black with the best fit in red. The continuum (gray dashed line) 
and individual emission-line components (grey gaussian curves) are also in- 
dicated respectively. In lower panels, the residuals are shown with same scale 
as that of upper panel. 

determination of the continuum characterized by a power law, 
fx cc \~ a . 

We employ multiple Gaussian components to describe the 
line profile. It is common practice to make such an assump- 
tion on the intrinsic shape of individual components, even 
though it has been demonstrated that broad-emis sion fines in 
AGN are not necessarily of such a shape (e.g., ICollin et"aT1 
2006). The Ha line is fit with two or three Gaussians (includ- 
ing a narrow component) and the [N II] A6548,6684 lines with 
a pair of Gaussians. The ratio of the [N II] lines is fixed at the 
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FIG. 4. — Examples of our fits for the Mg II line for the same AGN as 
shown in Figure|3] The top panel shows the spectral range around Mg II. The 
best fit model is indicated as a red solid line. The different components are 
as indicated: Fe-emission (blue), pseudo continuum (purple), and Gaussian 
components (green and brown). The residual is shown in the lower panel. In 
the upper right corner for each panel, we show the black hole mass distribu- 
tion computed from our Monte Carlo tests based on the uncertainties of the 
line width and continuum luminosity measurements. 



laboratory value of 2.96. The narrow width of the [N II] lines 
is fixed to match the narrow component of Ha. The width of 
the narrow components is limited to 420-800 km s" 1 (a range 
not corrected for intrinsic dispersion). The velocity profile of 
the broad components is characterized by the FWHM, mea- 
sured using either one gaussian or the sum of two gaussians. 
We then correct the velocity width for the effect of instrumen- 
tal dispersion to achieve an intrinsic profile width. The Ha 
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luminosity discussed throughout this work is the sum of the 
broad components. 

There are cases for which the fitting routine returns a so- 
lution with the width of the narrow component pegged at the 
upper bound of 800 km s" 1 . It is worth highlighting that this 
minimization routine stops the fitting procedure when a pa- 
rameters hit a limit thus the returned values are not the true 
best-fit values. For these, we inspect all fits by eye and de- 
cide whether such an additional broad component is real. For 
many cases, we can use the [O III]A5007 line profile, within 
the FMOS spectral window, to determine whether such a fit 
to the narrow line complex is accurate. In addition, we can 
use the available optical spectra for such comparisons. When 
the level of significance of the narrow line is negligible, we 
rerun the fitting routine and fix the narrow line width to the 
spectral resolution of FMOS, ~ 420 km s" 1 . In Figure [3] we 
show three examples of our fits to the Ha emission line that 
span a range of line properties. 

3.2. Mg II 

We fit the Mg I I em ission line, as done in 
ISchramm & Silvermanl (1201 2b . o bserved in optica l spec- 
tra primarily from zCOSMOS j LiUy_etaL| |20 Q9|), Mag- 
ellan/IMACS dTrump et alJ 120071) . and SPSS (lYorketal.1 
2000). The emission line is modeled by a combination of 
one or two broad Gaussian functions to best characterize 
the line shape. We first remove the continuum (before 
attempting to deal with the emission lines) by fitting the 
emission in a window surrounding the line. As with Ha, 
a power-law function is chosen to best characterize the 
featureless, non-stellar light attributed to an accretion disk. 
We further include a broa d Fe emission component b ased on 
an empirical template (Vestergaard & Wilkes] 12001 1) that is 
convolved by a Gaussian of variable width and straddles the 
base of the Mg II emission line. A least square minimization 
is implemented to determine the best-fit parameters. When 
possible, we optimize residuals of the fits on a case-by-case 
basis by trying to minimize the number of components. 
Absorption features are either masked out or interpolated 
across. The fit returns two parameters required for black hole 
mass estimates: FWHM and monochromatic luminosity at 
3000 A. Examples of our fits to the Mg II line are presented 
in Figure |4] 

4. RESULTS 

We can determine how closely the parameters (i.e., lumi- 
nosity and FWHM) required to estimate single-epoch black 
hole masses agree between the Ha and Mg II emission lines. 
Any systematic offset or inherent scatter may only add ad- 
ditional uncertainty to the derived masses. We essentially 
want to establish whether or not the kinematics of the BLR 
is consistent with photoionized gas in virial motion around 
the SMBH. We provide all measurements and derived masses 
in TableQ] 

4.1. Luminosities 

Our first concern is to determine whether the Ha emission- 
line luminosity scales appropriately with the UV continuum 
luminosity. For the following analysis, we do not correct 
for extinction due to dust and any contamination by the host 
galaxy; the impact of these, thought to be small, will be quan- 
tified in a later study. In Figure [5] the continuum luminos- 
ity, XL\ at 3000 A, is plotted against the emission-line lumi- 
nosity of Ha. Our data (as shown by the red points) spans 




log L Ra [erg/s] 

FIG. 5. — Comparison between the monochromatic luminosity at 3000 
A and the Ha emission-line luminosity. AGNs in our sample are shown as 
red circles. The open circles mark published observations of SDSS quasars 
(grey circles: 0.35 < z < 0.40; black circles: 1.5 < z < 2.2). A linear fit to 
both our data only (red line) and the joint sample (black line) is given. The 
distribution of log(AL3ooo/£Ha) as a histogram is shown in the sub panel. 
The colors of their histograms are same as the symbols. Errors represent a 
la level of uncertainty. 

two decades in luminosity and exhibits a clear correspon- 
dence between continuum and line emission. Based on our 
AGN sample, we determine the best-fit linear relation to be 
log(AL 3 ooo) = (0.82 ± 0.08) x logL Ha + (9.3 1 ± 3.47). For the 
linear fitting, we adopt a FITEXY method dPress et al.ll 19921: 
iPark et al.l2012l) . The level of dispersion of the data about this 
fit is 0.20 dex that will contribute to the dispersion in the final 
mass estimates. 

We can compare our data set to more luminous quasars 
from the SDSS. In partic ular, we i d entify 327 quasars 
from the SDSS sample in IShen et alJ (1201 lb with 0.35 < 
z < 0.40, that cover a similar luminosity range as our 
high redshift sample. These were selected from 1178 
quasars at 0.35 < z < 0.40 based on high-quality data 
determined by adopting the following criteria: error in 
log(M BH /M ) < 0.5, log(FWHMHa/km s -1 ) > 2.9, and 
log(FWHMMgii/km s" 1 ) > 2.9. In additio n, we include data 
from a recent study by IShen & Qui (1201 2b that provides the 
emission line properties including Ha line of high-luminosity 
quasars from the SDSS with 1.5 < z < 2.2. These samples 
are added to our data shown Figure [5] We clearly see that 
SDSS quasars fall along the AL3ooo-£hq relation as estab- 
lished above. Furthermore, the SDSS quasars have similar 
dispersion at both low-z and high-z samples to our sample, 
a = 0.20 and a = 0.15, respectively. We highlight that our 
AGN sample nicely extends such comparisons between con- 
tinuum luminosity and line emission at higher redshifts and 
to lower luminosities. We are able to effectively establish a 
wider dynamic range, not present in the high-z SDSS sam- 
ple due to the limited luminosity range around log(AZ/;ooo) ~ 
46.2. By merging all three samples, we find the following 
relation based on a linear fit: log(AL3ooo) = (0.96 ± 0.01) x 
logL Ha + (3.00 ±0.53). 

While there is very good agreement between the UV con- 
tinuum and emission line luminosity, there is a small dif- 
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TABLE 1 

Catalog of emission-line properties and black hole masses 



CID a 


Field 


RA [J2000] 


DEC [J2000] 




log L [ei 


g s -1 ] 


log FWHM [kms" 1 ] 


logM BH [M ] 


Ha 


3000 A 


Ha 


Mg II 


Ha 


Mg II 


178 


COSMOS 


149.58521 


+2.05114 


1.350 


43.76±0.03 


45.35±0.03 


3.685±0.013 


3.811±0.079 


8.68±0.03 


8.78±0.16 


5275 


COSMOS 


149.59021 


+2.77450 


1.400 


44.10±0.13 


45.50±0.02 


3.835±0.028 


3.973±0.072 


9.18±0.09 


9.18±0.14 


322 


COSMOS 


149.62421 


+2.18067 


1.190 


43.16±0.14 


45.02±0.02 


3.440±0.044 


3.581±0.081 


7.84±0.12 


8.17±0.16 


192 


COSMOS 


149.66358 


+2.08522 


1.220 


43.24±0.14 


45.05±0.02 


3.339±0.075 


3.490±0.035 


7.68±0.17 


8.00±0.07 


157 


COSMOS 


149.67512 


+1.98275 


1.330 


42.98±0.13 


44.84±0.03 


3.701±0.040 


3.632±0.017 


8.28±0.11 


8.19±0.04 



3 The Chandra obs. ID. 



ference that may impact, even slightly, our comparison of 
the masses. Based on the FMOS sample, the mean ra- 
tio (log(AL 3 ooo/^Ho!)) is 1.54 ±0.03, slightly higher than 
that found for both the low-z and high-z SDSS quasar sam- 
ples mentioned above; (logfALsooo/^Ha)) = 1-40 ±0.01 and 
(log(AL3ooo / Ltia)) - 1-36 ±0.02, respectively. This can be 
seen in the inset histogram in Figure [5] If this was due to the 
effect of dust extinction, one would find the opposite trend 
with reduced UV continuum relative to the Ha line emission. 
There may be other explanations such as an unde rlying SEP 
that may be different for X-ray selected samples (El vis et all 
l2012tlHao et al.ll2~012h and has an impact on the response seen 
in photoionized gas, or the effect of the host galaxy on the 
aperture corrections that differs in each band. While this is- 
sue is of importance, we reserve a detailed investigation to 
subsequent work since it is beyond the scope of this paper to 
adequately demonstrate such effects. Here, we are primar- 
ily concerned with the magnitude of a luminosity offset and 
whether it contributes to an offset between the masses. With 
black hole mass scaling with the square root of the luminos- 
ity (see below), the offset in luminosity, as determined above, 
amounts to a very small offset in log(MBH/A^o) of 0.07. 

4.2. Velocity widths 

A second pillar for the use of Mg II as a black hole mass 
indicator is that the emitting-line gas is located essentially 
within the same clouds that emit Balmer emission. While 
there are claims that this is the cas e by comparing the ve- 
locity profile of M g II with H/3 (e.g.. lMcLure & Jarvisll2002t 
IShen et all 1201 lb , there are r eported differences and trends 
that are not well understood (IWang et al.l 120091 IShen & Liul 
l2012tlTrakhtenbrot & Netzerl2012l) . For example, Mg II tends 
to be narrower than H/3 with a differenc e significantly larg er 
at higher velocity widths (see Figure 2 of lWang et al.1 12009). 

Our aim here is to compare the FWHM of the Mg II and Ha 
emission lines using a sample not yet explored, namely the 
moderate-luminosity AGNs at high-redshifts in survey fields 
such as COSMOS. In Figure [6] we plot the emission-line ve- 
locity width between Ha and Mg II emission lines. Based on 
the FMOS sample only, a positive linear correlation is seen 
between the velocity width of the two emission lines with 
the mean ratio of (log(FWHM M gii/FWHM H a)) = 0.089 ± 
0.022 (a = 0.143). Our data is in very good agreement with 
those from the SDSS, i.e., (log(FWHM Mg ii/F WHMHa)) = 
0.033 (a = 0.176) for the low-z sample from IShen et aTl 
d20Tl . (log(FWHMivi p it /FWHMHn)) = .00 (a = 0.076) for 
the high-z sample from S hen & Liul d2012l) and recent FMOS 
results from SXDS (see Nobuta et al. 2012). Based on a lin- 
ear fit to the data shown in Figure [6] we measure a slope that 
is consistent with unity; 0.898 ±0.132 for FMOS only and 
1.001 ±0.001 for FMOS+SDSS, and that cannot substantiate 
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FIG. 6. — Comparison between the FWHM of Ha with that of Mg II. 
Red, gray, and black circles are the same as those in Figure [5] A linear fit to 
both our data only (red line) and the joint sample (black line) is very similar 
to a one-to-one relation indicated by the dashed line. The distribution of 
log(FWHMMgii/FWHMHa) is also shown in the sub panel. 



the claim by lWang et al.l (120091) for a shallower value. These 
results are supportive of a scenario where the Mg II and the 
Ha emitting regions are essentially co-spatial with respect to 
the central ionizing source. 

There are a few noticeable outliers well outside the disper- 
sion of the sample. These objects then appear as outliers when 
comparing their masses based on different lines in the next 
section. Upon inspection, we find that these are the result of 
the FMOS spectra having low S/N. In some cases, there may 
be a fit to the Ha emission line, based on different parameter 
constraints, that is equally acceptable to the original fit as as- 
sessed by a chi-square goodness of fit and has a velocity width 
in closer agreement with Mg II. Although, we refrain from 
such selective fitting in order to present results that may be 
obtained from using similar fitting algorithms on larger data 
sets where such close inspection is not feasible. 

4.3. Virial single-epoch masses 

We now calculate black hole masses (Mbh) based on our 
single-epoch spectra using both (i) L3000 and FWHMMg 11, and 
(ii) Lua and FWHMh q . This calculation can be expressed as 
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o Shen & Liu (2012) 




log ^ Ha [M Q ] 



FIG. 7. — The comparison of the black hole mass estimated by Ha line 
with that by Mg II line. Red, gray, and black circles are the same as those in 
Figure[5] The fit to both our data (red line) and the joint sample (black line) 
is nearly equivalent to a one-to-one relation (dashed line). The distribution of 
log(MMgii/MH a ) is also shown in the sub panel. 



follows: 

, ,'Mbh 
log 



/ XL X orLune 
5 V10 44 ergs s" 1 



M Q J "°V!0 44 ergss-7 "°Vkm: 

We ex plic itly use the recipes pro vided by iGreene & Hoi 



/ FWHM \ 
+ clog( - 

111 



We ex plic itly use trie recipes pro vided by Ureen e & riol 
(2005) and McLure & Dunlop] (12004J) for the cases of Ha and 
Mg II lines, respectively: 

(a, b, c) = (1 .221, 0.550, 2.060) for Ha, 
(a, b, c) = (0.505 , 0.620, 2.000) for Mg II. 

We note that the calibration of the rel ation for Mg II has been 
carried out by many st udies (e.g., fVestergaard & Petersonl 
2006: lMcGilletaiT2 008) and there are known differences be- 
tween them. 

In Figure [7] we show Mbh for our sample derived 
from Ha and Mg II. Our sample spans a range of 7.2 < 
log(MBH/A^o) ^9.5 consistent with that rep orted by pre- 
vious studies of type 1 AGNs in COSMOS dMerloni et alj 
l2010tlTrump et al.l201 ll) and is complementary to the higher- 
L quasar sample at similar redshifts with log(MBH/M ) > 9 
dShen & Liull20l2b . We find the average (dispersion) in the 
black hole mass ratio of (log(Myi g n/Mn a )) is 0.17 (a = 0.32) 
for our FMOS sample. These results are similar to that deter- 
mined from the SDSS sample; the average (dispersion) in the 
black hole mass ratio is -0.05 (a = 0.39) at low z and -0.03 
(a = 0.20) at high z. While an offset of 0. 17 dex is seen in the 
FMOS sample, we conclude that the recipes established us- 
ing local relations give consistent results between Mg II- and 
Ha-based estimates. 



5. CONCLUSIONS 

We have investigated the emission line properties of AGNs 
in COSMOS and ECDF-S to establish whether Mg II and 
Ha provide comparable estimates of their black hole mass. 
This study is the first attempt to do so for AGNs of moderate 
luminosity, hence lower black hole mass (10 7 M Q < Mbh < 
1O 9 M ), at high redshift that complements studies of more 
luminous quasars. Our results clearly show that the velocity 
profiles of Mg II and Ha are very similar when characterized 
by FWHM and the relation between continuum luminosity 
and line luminosity is tight. We then find that virial black hole 
masses based on Mg II and Ha have very similar values and a 
level of dispersion (a ~ 0.3) comparable to luminous quasars 
from SDSS. It is important to kee p in mind that these re- 
sults pertain to spec ific calibrations (Mc Lure & Dunlopl2004t 
IGreene &~H o 2005) for estimating black hole mass. Th e use 
of other recipes, such as provided by tMcGill et al.ld2.Q08b . will 
show a discrepancy larger than seen here. To conclude, the 
locally-calibrated recipes for black holes masses using Mg II 
and Ha are applicable for fainter AGN samples at high red- 
shift. These results further support a lack of evolution in 
the physical properties of th e broad line region in terms of 
quant ities such as a ox (e.g., Steffen et al. 20061 : iGreen et ail 
120091) . emission-lin e strengths (e.g.. lFanll2006l). and the in- 
f erred metallici ties (iMaiolino et al 1 120011: iNagao et a"Dl2006t 
iMatsuokaet al.ll2011l) . 

As a final word of caution, such estimates of black hole 
mass are likely to have inherent dispersion as discussed above 
and systematic uncertainties that are not yet well understood. 
For instance, recipes for estimating black hole mass depend 
on the assumption that the gas is purely in virial motion. This 
is unlikely to be true for all cases since both outflows and in- 
flows are common in AGNs. Even so, there is evidence that 
the virial product of mass and luminosity (as a proxy for the 
radius to the BLR) is a useful probe of the central gravita- 
tional potential. In the very least, it is important to establish 
the level of dispersion in such relations since observed trends 
usually rely on offsets comparable to the dispersion such as 
the redshift evolution o f the relation between black holes and 
their host galaxies (e.g. JPeng et al.ll2006t lMerlonietaT1l2010t 
ISchramm & Silv erman 20121). 
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